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Abstract The angular distribution of B — KU (£ = e, ji, x) 
depends on two parameters, the lepton forward-backward 
asymmetry, Apg, and the flat term, F^. Both are strongly 
suppressed in the standard model and constitute sensitive 
probes of tensor and scalar contributions. We use the lat¬ 
est experimental results for £ = ji in combination with the 
branching ratio of > /i/r to derive the strongest model- 
independent bounds on tensor and scalar effective couplings 
to date. The measurement of provides a complementary 
constraint to that of the branching ratio of B^ —> /i/r and al¬ 
lows us—for the first time—to constrain all complex-valued 
(pseudo-)scalar couplings and their chirality-flipped coun¬ 
terparts in one fit. Based on Bayesian fits of various sce¬ 
narios, we find that our bounds even become tighter when 
vector couplings are allowed to deviate from the standard 
model and that specific combinations of angular observables 
in B —> A* are still allowed to be up to two orders of mag¬ 
nitude larger than in the standard model, which would place 
them in the region of LHCb’s sensitivity. 


1 Introduction 

With the analysis of the data collected by the LHCb Collab¬ 
oration during run I at the Large Hadron Collider (LHC), we 
now have access to rather large samples of rare B-meson de¬ 
cays with branching ratios below 10^^. As a consequence, 
angular analyses of three- and four-body final states can 
be used to measure a larger number of observables than 
previously possible at the B factories BaBar and Belle. In 
this work we focus on rare B decays driven at the parton 
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level by the flavor-changing neutral-current (FCNC) transi¬ 
tion b —> s££ that constitutes a valuable probe of the standard 
model (SM) and provides constraints on its extensions. 

The angular distribution of B —^ K ££—normalized to the 
width f]—in the angle 9e between B and £^ as measured in 
the dilepton rest frame is 

LHCb analyzed their full run 1 data set of 3 fb^' and mea¬ 
sured the angular distribution of the mode B+ — K+fi^, i.e. 
£ = Id [1], with unprecedented precision. They provide the 
lepton-forward-backward asymmetry Apg and the flat term 

in CP-averaged form and integrated over several bins in 
the dilepton invariant mass q^. Similarly, the CP-averaged 
branching ratios, = PbI]i, [2] and the rate CP asymme- 
try A^p [3] are also available from 3 fb^'. 

Both angular observables, F^ and Apg, exhibit strong 
suppression factors for vector and dipole couplings present 
in the SM, thereby enhancing their sensitivity to tensor and 
scalar couplings [4, 5]. A similar enhancement of scalar cou¬ 
plings compared to helicity-suppressed vector couplings of 
the SM is well-known from B^ —^ /i/i. Unfortunately the 
limited data set of B K*{^ K7t)££ from LHCb [6] did 
not yet allow to perform a full angular analysis without the 
assumption of vanishing scalar and tensor couplings in this 
decay mode. In the future with more data or special-purpose 
analysis techniques like the method of moments [7], certain 
angular observables in B — K*££ will provide additional 
constraints on such couplings, as for example J(,c [8] and 
the linear combinations {J\s — iJis) and {Jxc + Jic) [5, 9] as 

( 2 ) ( 3 ) 

well as the experimental test of the relations FIj = FIj and 
Jj =0 [5] at low hadronic recoil. 

Here we exploit current data from B+ —Refill and 
Bi —> /i/r to derive stronger constraints than before on tensor 
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and scalar couplings in various model-independent scenar¬ 
ios and study their impact on the not-yet-measured sensi¬ 
tive observables in Z? — K*ii. In Section 2, we specify the 
effective theory of \AB\ — 1 decays on which our model- 
independent fits are based. Within this theory, we discuss the 
dependence of observables in B —> KU and B —> K*U on the 
tensor and scalar couplings in Section 3 and specify also the 
experimental input used in the fits. The constraints on tensor 
and scalar couplings from the data are presented for several 
model-independent scenarios in Section 4. Technical details 
of the angular observables in B —?> K*££, the branching frac¬ 
tion of Bs —> / 2 /r, the treatment of theory uncertainties, and 
the Monte Carlo methods used are relegated to appendices. 


2 Effective Theory 


In the framework of the \AB\ = |ziB| = 1 effective theory 
Cj (flfo) 1^7 + Cji {Bb)^T 
-fh.c., 


4Gf ae * 


+ L 


( 2 . 1 ) 


the most general dimension-six flavor-changing operators 

(i) 

’ mediating b ^ sy and b —>■ s££ are classihed accord¬ 
ing to their chiral structure. There are dipole {i = 7,7') and 
vector {i = 9,9', 10,10') operators 

^9{9')=[^ynPL(R)b][lf^], ( 2 . 2 ) 

^fo(lO') = [s7tiPL(R)b][£ty5£] , 
further scalar {i = S,S\P,P') operators 

Ks')=[^PRiL)b]m, 

(2.3) 

and tensor (i = T,T5) operators 
= [sa^yb] [£a>^''£] , 

(2.4) 

^1-5 = [sa^yb] [f(T^''75f] , 

where the notation [sayiv^] is also 

used frequently in the literature. The respective short- 
distance couplings, the Wilson coefficients are 

evaluated at a scale of the order of the h-quark mass fib ~ mb 
and can be modified from SM predictions in the presence of 
new physics. 

The SM values '^7910 are obtained at next-to-next-to 
leading order (NNLO) [10, 11] and depend on the funda¬ 
mental parameters of the top-quark and VT-boson masses, 
as well as on the sine of the weak mixing angle. Moreover, 


they are universal for the three lepton flavors £ = e, fl, 1. 
All other Wilson coefficients are numerically suppressed 
or zero; — ms/mb^f^, ^ mbini/m^, and 

”^ 9 ' TT 5 “ Wilson coefficients of the four-quark 

current-current and QCD-penguin operators as well as of the 
chromomagnetic dipole operators are set to their NNLO SM 
values at fib =4.2GeV [10, 11]. 

For the rest of this article, we will suppress the lepton- 
flavor index on the Wilson coefficients —7 and opera¬ 
tors ^ In Section 4 we exploit data with £ = fi only, 
hence all derived constraints apply in principle only to the 
muonic case but can be carried over to the other lepton fla¬ 
vors £ = e,T for NP models that do not violate lepton flavor. 
In general, the Wilson coefficients are decomposed into SM 
and NP contributions -f but often we will 

use ^ for Wilson coefficients with zero (or suppressed) SM 
contributions synonymously with 


3 Observables and experimental input 

The full dependence of and Apg on tensor and scalar 
couplings has been presented in [4, 5], adopting the effec¬ 
tive theory ( 2 . 1 ), i.e. neglecting higher-dimensional opera¬ 
tors with dim > 8 . These results imply that for SM values of 
the effective couplings 

4b(?")|sm = 0- (3.1) 

Hence, for £ = e, fl both observables are quasi-null tests. 
The flat termB)^(g'^)|sM is strongly suppressed by small lep¬ 
ton masses for the considered kinematic region I <q^< 22 
GeV^ [4, 12, 13]. Nonzero values of ApgjsM can be in¬ 
duced by higher-order QED corrections, which will mod¬ 
ify the simple cos 9i dependence of the angular distribu¬ 
tion ( 1 . 1 ), however, currently there is no solid estimate avail¬ 
able for this source of SM background*. This picture does 
not change in the presence of new physics contributions 
to vector and dipole operators / = 7,7',9,9', 10,10'. On the 
other hand, nonvanishing tensor or scalar contributions are 
enhanced unless the dynamics of NP implies similar sup¬ 
pression factors, i.e., lepton-Yukawa couplings for or Ue 
in the case of Apg. In particular, F^ is very sensitive to tensor 
couplings (see (3.3) below) and Apg is sensitive to the inter¬ 
ference of tensor and scalar couplings (see (3.6) below). 

* Logai'ithmically enhanced NLO QED correetions to B —> Xs££ [14] 
turn out to be non-negligible for angular observables, however, anal¬ 
ogous corrections to 6 —> KU are partially included in the event sim¬ 
ulation of the experimental analysis [1] via PHOTOS [15]. As pro¬ 
posed in [16], the measurement of higher moments of the decay dis¬ 
tribution (1.1) could give an estimate of the size of higher-order QED 
corrections, but still admixed with contributions of dim > 8 operators. 
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Fig. 1 The sensitivity to the tensor coupling Re('^j.) of and SSfi in —> K^p.^i as well as {Jic +J 2 c), and {Ju — 3J2.<) in 

Angular observables are rescaled by the lifetime of the B meson, Tg. The bands represent the theory uncertainties at 68% and 95% probability of 
the prior predictive. Two sets of bands are shown for = 0 (blue) and = —1.1 (red). If available, the gray band indicates the latest 68% 
confidence interval reported by LHCb. All observables are integrated over 6 bins denoted as (...) j ^2 ^ 2 ] to match LHCb. 
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There are some angular observables 7, in B > /T* (—)■ 
Kn)M with the same properties; i.e. tensor and scalar con¬ 
tributions are kinematically enhanced by a factor \f(^jmi 
over vector ones present in the SM or their respective in¬ 
terference terms. These are J(,c and the two linear combi¬ 
nations (7is — 372i) and (Jic+Jic) with explicit formulas 
given in Appendix A. In our fits and predictions we include 
all kinematically suppressed terms. But for the purpose of 
illustration, we now consider the analytical dependence for 
vanishing lepton mass. In this limit, 

Jec Re[(<4 - ^p,)^T - (3-2) 

is sensitive to the interference of tensor and scalar operators, 
complementary to Apg in B ^ KU 


A^b oc Re[(<4 -f -f -f In, (3.3) 

i.e., with an interchange of tensor contributions T o 75. We 
note also that J^c contributes to the lepton forward-backward 
asymmetry of B —>■ K*££ being {J(,s+J(,c/2). Since it has 
to compete with in this observable, a separate measure¬ 
ment of J^s and J(,c is necessary. 

Only tensor contributions enter 

(/ 1 .- 372 .)- + . ( 3 - 4 ) 


where the dots indicate different kinematic and form-factor 
dependencies. But tensor and scalar contributions enter 


(7ic -I-72 c) + I'^rsP) 

+ • • • ( 1^5 ~ ) 


(3.5) 


which is similar to the dependence of in B — K££ 




... + 

-f . . . ( I -f '^^1 I ^ + \^ p ~\- '^ p ! I 


(3.6) 


'Tf. 


Concerning the involved kinematic factors—see [4, 
5] —are such that tensor and scalar couplings contribute 
only constructively/cumulatively, apart from cancellations 
among and ‘^gi^^piy Interference terms in the numer¬ 
ator of of the form and i^pp, x jg,) 

are suppressed by mij They become numerically rel¬ 
evant in case ^ ^ or "^pp, '^jg jg/ where the 

smallness of pp, is of the same level as the suppression 

factor mil accompanying the large vectorial SM Wil¬ 
son coefficients ±4. This implies, however, no large 

enhancement of F^ over the SM prediction. 

On the one hand, the observables F^ (3.6) and Apg (3.3) 
are measured in the angular distribution (1.1) of B —> K££ 
normalized to the decay width T^ such that uncertainties 
due to form factors can cancel in part [4, 5]. On the other 


hand, J(,c, [Ju — and (7ic + 72e) appear in the unnor¬ 

malized angular distribution of B —> K*££. “Optimized” ver¬ 
sions Bi, Ml, and M 2 for the low-g'^ region for which form 
factors cancel in the limit of mh °° have been identified 
in [9]. For the high-q'^ region, potential normalizations are 
discussed in Appendix A, which could serve to form opti¬ 
mized observables for special scenarios of either vanishing 
chirality-flipped vector or tensor or scalar couplings. In the 
most general case, however, there are no optimized observ¬ 
ables at high q^. Although form factors do not cancel in this 
case, it might still be preferable to use normalizations, for 
example when the overall normalization of B ^ form 
factors constitutes a major theoretical uncertainty. 

To illustrate the sensitivity of F^ to tensor couplings, we 
compare it in Figure 1 to the branching ratios of B —> *)££ 
for £ = jJ., integrated over one low-q^ and one high-q'^ bin. 
The details of the numerical input and the uncertainty propa¬ 
gation can be found in Appendix B and Appendix C. In light 
of the hint of new physics in from recent global analy¬ 
ses of b s{y, ££) data [26-30], we show predictions for 
= —1.1 in addition to = 0. 

From Figure 1, the highest sensitivity to tensor couplings 
of any B K observable is attained by at high q^ due 
to a partial cancellation of form factors [5]. If the experi¬ 
mental uncertainty could be reduced further, B^ would give 
a very strong constraint on a simultaneous negative shift in 
and The prediction of —> KfLjj.) is essentially 

insensitive to but sensitive to A stronger impact on 
global fits, however, would require a reduced theory uncer¬ 
tainty. 

The observable £^{B K*p.iJ.) shows moderate depen¬ 
dence on at least at low q^ and has some impact on the 
constraints on tensor couplings as will be discussed in Sec¬ 
tion 4. At the moment, theory and experimental uncertainty 
are of similar size. 

{Jic + Jic) is sensitive to in both q^ regimes. At 

low q^, it is mildly affected by , whereas at high q^ it is 

unaffected. Regarding {J\s — 3J2s), the situation is reversed: 
here the strong dependence on ‘^p appears at low q^. Over¬ 
all, F^, {Jic + Jic) at high q^, and (7ij — 3J2s) are sensitive 
to 'rfp and theoretically very clean around ‘^p = 0. 

From the available measurements, F^ at high q^ cur¬ 
rently provides the most stringent constraints on the size of 
tensor couplings. Moreover, the dependence on vector cou¬ 
plings is such that < 0 leads to stronger constraints 
on ^p than « 0. 

Important additional constraints on scalar couplings 
come from the branching ratio of B^ —> (2/4 as given in 
(A.8). It provides the most stringent constraints on the mod¬ 
uli I and — 'Wp,\ and further depends only on 

("^10 ~ '^ 10 ')- Thus it is complementary to B^ in B —^ Kit, 
see Eq. (3.6). 
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Channel 

Constraints 

Kinematics 

Source 


B,i —> pp 

^ = /dr^(T) 

- 

[17-19] 



% 

6 [1,6], [14.18, 16], [> 16]GeV2 
?2 6[l.l,6.0],[15.0,22.0]GeV2 

[20] 

[2] 

B+ -i-K+pp 

^FB 

p2 6[l.l,6.0],[15.0,22.0]GeV2 

[1, 20] 



6 [1.1, 6.0], [15.0, 22.0] GeV^ 

[1] 


^CP 

6 [1.1, 6.0], [15.0, 22.0] GeV^ 

[3] 


% 

6 [1,6], [14.18, 16], [> 16] GeV^ 

[20-22] 

B®-)-K*°pp 


6 [1,6], [14.18, 16], [> 16] GeV^ 

[20-22] 


^CP 

6 [1.1, 6.0], [15.0, 22.0] GeV^ 

[3] 

B —> K form factors 

fo.+.T 

= 17, 20, 23 GeV2 

[13] 

B —> K* form factors 

L, Ao,1,2, 7’p2,3 

L, Ao, 1,2, 7i.2,3 

1/2= 0.1,4.1, 8.1, 12.1 GeV^ 

6 [11.9, 17.8] GeV^ 

[23] 

[24, 25] 


Table 1 List of all observables of the various b ^ sp.^1 decays entering the fits with the respective kinematics and experiments that provide the 
measurements. LCSR and lattice results of S —> K^*'l form factors are used to constrain a i^^-dependent form-factor parametrization. For more 
details see Section 3 and Appendix B. 


Eventually we also explore the effect of interference 
with NP contributions in the vector couplings g, jq/ 
on the bounds on tensor and scalar couplings. For this pur¬ 
pose we include also the branching ratio, the lepton forward- 
backward asymmetry, and the rate CP asymmetry of B —> 
K*iiii as they provide additional constraints on the real and 
imaginary parts of g, jg jg;. The experimental input of all 
observables entering our fits is listed in Table 1 together with 
input for the B —!• form factors. More details on the lat¬ 

ter can be found in Appendix B. 


4 Fits and constraints 

There are no discrepancies between the latest measurements 
for £ = ji (throughout this section) of and Apg in B+ -> 
K^Bll and their tiny SM predictions; cf. Figure 1. Thus our 
main objective is to derive constraints on tensor and scalar 
couplings through the enhanced sensitivity of both observ¬ 
ables to these couplings compared to vector couplings. For 
this purpose, we will consider several model-independent 
scenarios, progressing from rather restricted to more gen¬ 
eral ones in order to asses the effect of cancellations due to 
interference of various contributions. 

For each coupling that we vary in a fit, we remain as 
general as possible, treat it as a complex number and use the 
Cartesian parametrization assuming uniform priors for ease 
of comparison with previous studies. Specihcally, we set 

Re('^^y pp/Pj-j) G [—1,1], 

Re('4,9Mo,io') e [-7,7], (4.1) 


and the same for the imaginary parts. The priors of the nui¬ 
sance parameters are given in Appendix B. 

We start with the scenario of only tensor couplings and 
see that they are well constrained by pj^ alone. In a sec¬ 
ond scenario we consider only scalar couplings in order to 
investigate the complementarity of P^ and > /i/i. Here 
we hnd that—for the first time—all complex-valued scalar 
couplings can be bounded simultaneously by the combina¬ 
tion of both measurements. Finally we consider as a special 
scenario the SM augmented by dimension six operators as 
an effective theory of new physics below some high scale 
Anp assumed much larger than the typical scale of elec- 
troweak symmetry breaking. In addition, the model con¬ 
tains one scalar doublet under SU{2)l as in the SM. For 
each scenario, we also investigate interference effects with 
new physics in vector couplings ‘^g g, jg jg,. Finally, we con¬ 
clude this section with posterior predictions—conditional on 
all experimental constraints—of the probable ranges of the 
not-yet-measured angular observables J(,c, (Jic+Jic) and 
{Jis — 3J2 s) in B K* fill and A^j- in ^ fifL. 

4.1 Tensor couplings 

In a scenario with only complex-valued tensor couplings 

J. 5 , the experimental measurement of pj^ constrains the 
combination + [‘^t’sP, up to some small interference 
of with vector couplings y, g g,; cf. (3.6). The ac¬ 
cording 68 % (95%) ID-marginalized probability intervals 
are listed in the second column of Table 2. From the third 
column, it is seen that the constraints become tighter when 
utilizing all observables in Table 1, mainly due to the sen- 
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data set 

set of couplings 
credibility level 

only Ff 

7'5 

68%, 95% 

Ff + other 

68%, 95% 

Fp + other 

75,9,10 

68%, 95% 

Re'^j. 

[-0.32,0.16], [-0.52,0.35] 

[-0.29, 0.09], [-0.43, 0.27] 

[-0.23,0.14], [-0.39,0.30] 


[-0.25,0.24], [-0.44,0.44] 

[-0.19,0.21], [-0.37,0.35] 

[-0.17,0.22], [-0.33,0.37] 

Re'^j.j 

[-0.25,0.24], [-0.44,0.44] 

[-0.19,0.17], [-0.33,0.33] 

[-0.18,0.16], [-0.32,0.32] 


[-0.24,0.25], [-0.44,0.45] 

[-0.21,0.19], [-0.37,0.35] 

[-0.20, 0.18], [-0.36,0.35] 

\Vt\, 

[0.13,0.43], [0.03,0.57] 

[0.09, 0.33], [0.02, 0.44] 

[0.10, 0.33], [0.02, 0.43] 


Table 2 The constraints on complex-valued when using measurements of only and other data in Table 1 (except > fJ-B)), 

and finally allowing complex-valued new physics contributions to jq. 


sitivity of the branching ratio of B —> K*jlii to tensor cou¬ 
plings (see also Figure 1). The latter stronger bounds are 
driven by the new lattice results of B ^ K* form factors that 
predict values above the measured ones [31]. Since tensor 
couplings contribute constructively to ^{B K* jliJ.), large 
values are better constrained. In this scenario with vanish¬ 
ing scalar couplings, current measurements of Apg(B+ — 
K^Bll) barely provide any constraint; cf (3.3). 

We also perform a fit with nonzero order to as¬ 

sess the robustness of the bounds with respect to interfer¬ 
ence. Note that -j, appears in linear combinations with 
g, such that its interference with tensor and scalar cou¬ 
plings is captured implicitly by allowing new physics in 
g,. Thus we fix y to the SM value without loss of gen¬ 
erality. In this case, by itself still provides bounds on 
I'^j. j-jl that are weakened by a factor of two since does 
not pose constraints on (in the chosen prior range). 

Once additional experimental measurements of Table 1 are 
taken into account, the potential destructive effects of new 
physics in jg become reduced and almost the same con¬ 
straints on are recovered, as shown in the last column 
in Table 2. If in addition we allow y pp, ^ 0 (not shown in 
Table 2), the credible regions further shrink by about 10%, 
which we attribute to the cumulative effect of pp, 7 ^ 0 
in fIj\ cf. (3.6). In summary, the measurement [1] of 
LHCb with 3 fb^* shrinks the previous bounds [5] on 
by roughly 50%. 

A keen observer may notice that in Table 2 the SM point 
= 0 is contained in every 68 % region in Cartesian co¬ 
ordinates but never in even the 95% region in polar coor¬ 
dinates^. This is a consequence of the general concentra¬ 
tion of measure. Another way to look at it is to transform 
the uniform prior density from Cartesian to polar coordi¬ 
nates. For the example of a single Wilson coefficient, say 
‘rfj, the transformed density is proportional to the determi¬ 
nant of the Jacobian which is |^p|. Since the Cartesian prior 
boundaries are much larger than the regions of high likeli- 

^But the bin with lower edge is always in the 99% region. 


hood, one could think the value of the boundary is irrelevant, 
but in fact it determines the peak of the prior in polar coor¬ 
dinates. In other words, the uniform prior on Re('^j,) and 
Im(^j.) favors larger values of Y^p\ even though we con¬ 
sider it consensus in the community that smaller rather than 
larger values are reasonable because = 0 in the SM. We 
suggest therefore that the default treatment be revised in the 
future to include available prior knowledge. 

4.2 Scalar couplings 

Scalar couplings y p p, enter F^ without kinematic 
suppression—see (3.6) —as the sum + whereas in 
the time-integrated branching ratio S§{Bs li) they ap¬ 
pear as the difference = S, P. Since the exist¬ 

ing measurement of Fj^ constrains the sum, the combination 
of and fill) allows us—for the first time—to 

bound the real and imaginary parts of all four couplings. The 
corresponding 2D-marginalized regions in the Re('^ ± If.,) 
{i — S,P) planes are shown in Figure 2. The corresponding 
plots for Im('^ ^re very similar to those shown and 

thus omitted. These bounds do not change when including 
all other data in Table 1, since the Apg(B — Kfift) requires 
interference of scalar with tensor couplings and other ob¬ 
servables are not very sensitive to scalar couplings. Quan¬ 
titatively, the constraint from ^lS{Bs —>■ ft ft) on — ^.,) is 
about a factor four to five stronger than the one of F^ on 
{% + %')■ 

Interference terms of ‘^p p, with vector couplings might 
weaken these bounds. For £§{Bs —t p-ll), the relevant term is 
(“^jg —(see (A. 7))andforF^itis (‘^ig + '^ig;) [4];both 
are suppressed by the factors nip/M b, and mpj respec¬ 
tively. Nevertheless, these terms become important for small 
'^p p, due to the large SM value of ~ —4.2. We compile 
bounds on complex-valued scalar couplings in Table 3 for 
only F^, only ^^{Bs P-B), their combination with all 

the other observables in Table 1 . Neither F^ nor P§{Bs li) 
alone can bound all four complex-valued scalar couplings. 








7 




Fig. 2 The constraints on complex-valued pp, from only (gray dotted), only 3§(Bs fill) (blue dashed), and the combination with all 

other data in Table 1 as well as nonzero Vg g/ jq jq/ (red solid) at 68% (darker) and 95% (lighter) probability. The constraints on Re('^p±'^p,) are 
identical to Re)”^^ ± apart from a small translation of the contours by (-1-0.2, +0.15). The SM predietion is indicated by the black diamond. 


data set 

set of couplings 
credibility level 

only Fp 

^S.S' .P.F 

68%, 95% 

only 3S{Bs llfl) 

68%, 95% 

all 

'^X,5',F.F'.9,9'.10.10' 

68%, 95% 

Re(‘rj-'ry) 

- 

[-0.10, 0.08], [-0.14,0.13] 

[-0.08,0.07], [-0.13,0.13] 


- 

[-0.07,0.07], [-0.11,0.12] 

[-0.07,0.07], [-0.12,0.11] 


[-0.36,0.39], [-0.69,0.68] 

- 

[-0.32,0.32], [-0.59,0.62] 

ImCWg + Vg,) 

[-0.37,0.35], [-0.68,0.66] 

- 

[-0.25,0.41], [-0.57,0.64] 

ReCTp-'Ty) 

- 

[0.05,0.20], [0.01,0.26] 

[0.00,0.16], [-0.07,0.22] 

ImC^p — ^p,) 

- 

[-0.07,0.08], [-0.12,0.12] 

[-0.07,0.09], [-0.14,0.16] 

R^C^P + Vp,) 

[-0.24,0.51], [-0.51,0.82] 

- 

[-0.12,0.52], [-0.42,0.78] 


[-0.36,0.37], [-0.67,0.67] 

- 

[-0.37,0.29], [-0.68,0.57] 


Table 3 The ID-marginalized constraints on complex-valued ^ pp, at 68% (95%) probability from measurements of only , only 3S{Bs —> 
and all the data in Table 1 and additional new physics contributions to Vg g, jq jq,. 


however their combination is capable to do so and more¬ 
over, the bounds are stable against destructive interference 
with vector couplings. In the case of jq jq, ^ 0, we find 
the following bounds with 68% (95%) probability; 

e [0.06,0.25] ([0.01,0.36]), 

|^5,|e [0.09,0.46] ([0.02,0.73]), 

\'^p\ e [0.05, 0.21] ([0.01, 0.32]), ' ’ 

\^p ,\G [0.11,0.42] ([0.02,0.63]). 

Allowing in addition ^ 0, the intervals of (4.2) are 
quite similar but in general (10-20)% narrower and shifted 
by that amount towards zero. Again, this can be explained 
by the cumulative effect of y pp, and in shown 
in (3.6). 

In the special case of real-valued couplings, ^(Bi —> iF) 
would lead to rings [32] instead of circles in Figure 2. Our 
results improve and extend previous bounds in the literature 


to the most general case of complex-valued couplings. For 
example they are a factor two to five more stringent than [33] 
and comparable to [34] once restricting to the simpler sce¬ 
narios considered there. 


4.3 SM-EFT-constrained scalar couplings 

In the following we consider a scenario in which it is as¬ 
sumed that there is a sizable hierarchy between the elec- 
troweak scale and the new-physics scale, Anp, and that the 
SM gauge symmetries SU{2)l xU{1)y are only broken at 
the electroweak scale. This results in the augmentation of the 
SM by dimension-six operators that respect the SM gauge 
group and are composed of SM fields only. Such a scenario 
becomes more and more viable for two reasons. The first is 
the discovery of a scalar resonance at the LHC in agreement 
with all requirements of the Higgs particle in the SM. The 



























second is the steadily rising lower bound on the mass of new 
particles reported by ATLAS and CMS in various more or 
less specific models. 

A nonredundant set of dimension-six operators of this 
effective theory (SM-EFT) that requires a linear realization 
of the electroweak symmetry was given in [35]. The match¬ 
ing of the SM-EFT to the effective theory of AB = 1 de¬ 
cays (2.1) at the scale jj. ^ mw of the order of the VT-boson 
mass was performed for vector couplings g jq in [36]. 
The matching of tensor (2.4) and scalar (2.3) operators [32] 
shows that SM gauge groups in conjunction with the linear 
representation impose the relations 


= = (4.3) 


on scalar couplings and require tensor couplings to be sup¬ 
pressed to the level of dimension-eight operators. In conse¬ 
quence only two scalar couplings ^ arise that scale as 
v/Anp)^ <C 1 where v niw denotes the scale of elec¬ 
troweak symmetry breaking. 

It must be noted that the relations (4.3) are a conse¬ 
quence of embedding the Higgs in a weak doublet along 
with the Goldstone bosons. For example, choosing a non¬ 
linear representation of the scalar sector allows additional 
dimension-six operators in the according effective theory, 
such that the couplings ^ pp, are all independent and ten¬ 
sor operators have nonvanishing couplings already at dimen¬ 
sion six [37]. 

Omitting for the sake of simplicity terms of order 
and OTj/ot*, the couplings y can be bound 

from [32] 

W(B,^U)oc + 

-— 2 —. 


In similar spirit, dropping terms of order m^jq^ and ms/mj, 
gives 




^^Re[(^,-<4,0 + ^,0')*] • 


(4.5) 


In the SM-EFT no relations between ‘^10 and'^jQ, arise, 
so they are in general additional independent parameters. 
Here we find that destructive interference with contribu¬ 
tions involving does not significantly alter the bounds 

on y. The results of two fits are shown in Figure 3. In 
the first ht, we set = 0 and include all constraints 

on B —i' Kfifi and —> /2/r. In the second fit, we allow 
*^ 10 %' 7^ ® further include all B K*fi^ constraints 
from Table 1. For both fits, all six 2D marginals of real and 
imaginary parts of vs. “^y have nearly circular contours 


of equal size that contain the SM point at the 68% level ex¬ 
cept for Re(‘^^) vs. Re(^y) where it is within the 95% cred¬ 
ible region. The regions hardly vary between the two fits. 

Since we consider here complex-valued couplings the 
allowed regions are circles rather than rings as for the case 
of real-valued couplings [32]. Compared to those rings, the 
circles are smaller because the probability moves from the 
ring towards the center of the circle. 


4.4 Tensor, scalar, and vector couplings 

The most general fit of complex-valued tensor and scalar 
couplings y pp, rp in combination with vector cou¬ 
plings g, JO jQ,—the combination of Section 4.1 and Sec¬ 
tion 4.2 —yields bounds very similar to those in tables 2 
and 3. The changes are only small and in fact the bounds 
tend to be even more stringent because tensor and scalar 
couplings can contribute to only constructively—see 
(3.6). As before, branching-ratio measurements of B —> 
K*p.ii help to improve the constraints on This 

demonstrates that even in the case of complex-valued cou¬ 
plings there is enough information in the data to bound all 
20 real and imaginary parts. 


4.5 Angular observables in B — K*£i 

Now we discuss what the fits tell us about likely values of 
observables that have not been measured yet but have sen¬ 
sitivity to tensor and scalar couplings. In the B —> A'*(— 
Kn)U decay, we again consider {J\s — 3J2s) and {Jic+J 2 c) 
as in Section 3 and additionally 7ge. We compute the poste¬ 
rior predictive distribution (see Appendix C) for each ob¬ 
servable integrated over the low-q^ bin [1.1,6] GeV^ and 
high-q'2 bin [15, 19] GeV^ matching LHCb’s range. The dis¬ 
tributions resemble Gaussians, thus we summarize them by 
their modes and smallest 68% intervals in Table 4 comparing 
the SM (prior predictive, = 0) to three NP scenarios. 
In each, we allow for interference with the vector couplings 
and additionally vary only (Section 4.1), only pp, 
(Section 4.2), and finally both tensor and scalar couplings 
(Section 4.4). 

We rescale 7, and combinations by the B'^-meson life 
time Tgo = 1.519 ps [38] to judge the experimental sensitiv¬ 
ity in the near future by comparing to current measurements 
of the branching ratio 

£^=^[2 (37i, - J2s) + (37i, - J2c)\ ■ (4.6) 

In the SM, the typical magnitude of the branching ratio of 
B —>■ K*i£ is approximately equal to 2 • 10 ^^ for both the 
€ [1-1,6] and [15, 19] GeV^ bins. For comparison, the 
predicted ranges for Tgo {Jis — 3 J 2 s) and Tgo (7ie + J 2 c) in the 
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Fig. 3 68% and 95% contours of the 2D-marginalized distributions of scalar couplings y in the scenario SM-EFT with all constraints in 
Table 1 (red) when marginalizing over nonzero Re, For comparison, we superimpose the corresponding contours from using only 

SS{Bs —> /ift) and all B —> Kfill constraints (dashed blue contours) with fixed ~ 0. The SM is indicated by the black diamond. 


observable 

ij^-bin [GeV^] 

SM 

T(5),9, 10 

S('),p('),9('),io(') 

r(5),9('),10('l 

TgO X 

[1.1,6] 

~0 

(0.6+1;*). lO-’’ 

(-0.1+2;*). 10-“> 

(0.2+};|).10-'' 

[15, 19] 

~0 

(2.1++>).10-'° 

(0.7^+^). 10-“ 

(-0.2+*;*). 10-'' 

TgO X Jlc) 

[1.1,6] 

(3.30+[];«).io-9 

(4.8+1;^). lO-** 

(2.6+°;*). 10-^ 

(2.8+1;*). 10-° 

[15, 19] 

(1.72+«;|6).10-'« 

(4.7+^). 10-^ 

(1.9^1];*). 10-'° 

(3.1^1*). 10-° 

%) X (Tls - 'iJ2s) 

[1.1,6] 

(2.44+0;46).10-'o 

(1.9+|;t).10-* 

(4.0+[;°).10-'° 

(1.5+1;?). 10-* 

[15, 19] 

(1.12+<];jS>).10-'0 

(7.1+5;^). 10-« 

(1.4^°;*). 10-'° 

(5.3^*;*). 10-° 


- 

1 

1 

[-1,1] 

[-1,1] 


Table 4 The posterior predictive 68% probability intervals of not-yet-measured angular observables for several new-physics scenarios given all 
the considered experimental constraints. The eorresponding values for the SM (prior predictive) are given, too, where 0” indicates zero in the 
considered approximation, see text for details. 


SM are suppressed by 2 — 3 orders of magnitude down to 
^(10-“*); cf. Table 4. 

The angular observable J(,c is strictly zero in the absence 
of tensor and scalar couplings. Nonzero contributions can 
be generated in the SM by QED corrections or potentially 
from higher-dimensional (d > 8) \AB\ — |zi5| = 1 operators, 
leading to parametric suppression by a^/ (Ttt) or nthing /. 
These factors should be compared to the potential suppres¬ 
sion present for tensor and scalar contributions in particular 
NP models in order to gauge their relevance. Our model- 
independent fits are still in a regime where such considera¬ 
tions are insignificant since current experimental measure¬ 
ments, in combination with theory uncertainties, do not yet 
impose sufficiently stringent constraints on tensor and scalar 
couplings. 

Beyond the SM, can become of order ^(10^®) 
in scenarios involving tensor couplings only and about 
^(10^***) in the presence of scalar couplings only. Both ef¬ 


fects are due to interference with vector couplings. Schemat¬ 
ically, Tgc is a function of x‘i^p +'if^ x niijx 

scalar x vector, and m^j x tensor x vector. The largest 

interval is obtained for the scenario without scalar couplings 
because then the uncertainty on the (tensor) couplings is 
largest. But even then, it seems that the experimental sen¬ 
sitivity will not be high enough to have an impact in global 
fits. 

Concerning {J\c + Jic) and — substantial de¬ 

viations from the SM prediction are again only possible in 
the presence of nonzero tensor couplings. In this case, an 
enhancement by two orders of magnitude is possible up 
to if?(10^^) at high (f- and also at low in the case of 
{Jis — We want to stress again that we make these 

statements conditional on all included experimental con¬ 
straints, the scenario, and our prior. In view of the current 
experimental precision of 20% on the branching ratio at 
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LHCb [21] with only 1 fb^\ corresponding to the ^(10^^), 
one can indeed hope for some sensitivity to such large ef¬ 
fects in {J\c + J 2 c) (Jis — ^Jis) for the not-yet-published 

3 fb^* data set. At least, we can hope for some measurement 
if the method of moments [7] is applied. 

For the CP asymmetry A^r{Bs Am) induced by the 
nonvanishing width of the Bs meson (cf. Appendix A), we 
find a rather uniform distribution in scenarios with nonzero 
scalar couplings. So any value in the range [—1, 1] is plausi¬ 
ble whereas the SM and the scenario with only tensor cou¬ 
plings predict a value of precisely one [39]; cf. the last row 
in Table 4. Hence any deviation from one would unambigu¬ 
ously hint at the presence of scalar operators. 


5 Conclusions 

We have derived the most stringent constraints to date on 
tensor and scalar couplings that mediate h —i/i/i transi¬ 
tions. They are based on the latest measurements of angu¬ 
lar observables and the lepton forward-backward asym¬ 
metry Apg in B+ from LHCb [1], supplemented 

by measurements of the branching ratios of Bg fiji and 

Both and Apg belong to a class of observables in 
which vector and dipole couplings—present in the stan¬ 
dard model (SM)—are suppressed (mostly kinematically by 
mil\f(^) with respect to tensor and scalar couplings. We 
provide predictions for the equivalent but not-yet-measured 
angular observables Tgc, {J\cPJic) and {Jis — 3J2s) in B 

In a Bayesian analysis of the complex-valued couplings 
of the effective theory, we find that the measurement of , 
especially at high-q'^, 

1 . imposes by itself constraints on tensor couplings 

such that the upper bound of the smallest 68% (95%) 
credibility interval is 0.43(0.57), superseding previous 
bounds. In combination with current data from B —> 
K*fxpL and lattice predictions of B —> F* form factors, 
the bounds are lowered to 0.33 (0.43), even in the pres¬ 
ence of nonstandard contributions in vector couplings. 

2 . for the hrst time allows to simultaneously bound all 
four scalar couplings y pp, due to its complementar¬ 
ity to ^{Bs —> /I/i). Even when taking into account de¬ 
structive interference with vector couplings, |^- -f | < 
0.3 (0.6) and -%,\< 0.1 (0.2) for i = S,P with at 
least 68% (95%) probability. Currently, the bounds from 
F^ are weaker than those from £§{Bs —/i/r) by about a 
factor of four. Future measurements of at LHCb and 
Belle II will further tighten the bounds. 

Moreover, measurements of (I = e) will provide con¬ 
straints on scalar couplings in the electron channel in the 


absence of a direct determination of the branching ratio 

of Bs ee. 

Our updated bounds on complex-valued tensor and scalar 
couplings are summarized in Table 2 and Table 3, account¬ 
ing also for interference effects with vector couplings. These 
bounds hold even in the most general scenario of complex¬ 
valued tensor, scalar, and vector couplings, showing that the 
data are good enough to bound the real and imaginary parts 
of all Wilson coefficients simultaneously. 

As a special case, we consider the scenario arising from 
the SM augmented by dimension-6 operators generalizing 
existing studies to the case of complex-valued couplings. In 
this scenario, tensor couplings are absent and additional re¬ 
lations between scalar couplings are enforced by the linear 
realization of the SU{2)l®U{\)y electroweak symmetry 
group. 

Our study of the yet unmeasured angular observables 
Jec, {J\c+J 2 c), and {Ju - 3J2s) in B ^ F*(-> Kk)bii (see 
Table 4) shows that despite the current bounds on tensor 
couplings, enhancements of up to two orders of magnitude 
over the SM predictions are allowed for {J\c + 72c) and 
(Jis — 3J2 s), placing them in reach of the LHCb analysis 
of the full run I data set. Our bounds on scalar couplings 
from Bs ^ p-B and F^, however, are already quite restric¬ 
tive permitting only small deviations from SM predictions 
in (7ic +72 c) and {Jis — 3J2s)- Notably, the CP asymmetry 
A^riBs pp) given nonzero scalar couplings can take on 
any value in the range [—1, Ij. 
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Appendix A; Angular observables in B —j K*U 

Here we focus on those angular observables 7, in B —> K*U 
in which tensor and scalar contributions are kinematically 
enhanced by a factor \f(^/mi compared to the vector con¬ 
tributions of the SM or the respective interference terms. The 
results and the notation follow [5]. 

For convenience, we split the time-like transversity am¬ 
plitude A, into two parts as 

Ai=Ai-\^Ap. (A.l) 

2 mi 
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Both terms depend on the scalar B ^ K* form factor 
Ao{q^), a normalization factor N, and the Kallen-function 
X (see [5]), but have a different dependence on the Wilson 
coefficients i — 10, 1Q\P,P': 


At=2N 


y/X 

y/¥ 


(^10 “^10'Mo, 


Ap = -2NVx f^~^^'] A,. 

(OTfe+OTj) 


(A.2) 


The lepton-flavor index i of Wilson coefficients is omitted 
for brevity throughout. 

In full generality, the angular observables in B —> K*i£ 
depend on seven transversity amplitudes with vector and 
dipole contributions, A^'^ || and A,, one scalar and one pseu¬ 
doscalar amplitude, As,p (^sp ~ ^s' p' ^^ tensor 

amplitudes, A||x,r±,o|| A,o,j||,ox ^he inter¬ 

esting combinations are 


-Zee =4)3^ Re 


2(A,oAj-A||^A*) 

+ [Mo +AoMJ -l-4A||j^A*] 


(A.3) 


resembles very much the branching ratio of the rare decay 
Bs —>■ M in the limit —>■ 







Mb ± nis 


-f 


'^p — ’^p ^ 2m^ 


nih + mj Mj 


+ JXT ('^10 “ ^lo) 


(A.7) 


Hence there is some similarity between {J\c +J 2 c) in B —> 
K*££ and ££§{Bs —^ ££) in their dependence on the couplings 
but the former has additional dependence on tensor and vec¬ 
tor couplings through the other transversity amplitudes. In 
B K*££, the helicity suppression factor A m^/q^ of vector 
couplings is weaker than the corresponding factor 4m^/ mJ 
in Bs —> ££. 

Concerning B^ —> ££, the expression (A.7) corresponds 
to the plain branching ratio at time f = 0. Due to the non¬ 
vanishing decay width ATJ, experiments measure the aver¬ 
age time-integrated branching ratio—denoted by ££§—and 
the two are related as [39] 


4 

3 


(7le + 72e) 


^(Ao+A^)+4A,o 


2 

+ 16j3MA|Ml" 


-f 



Ap 


2 

+ Pl\As\\ 


(A.4) 


3 ( 71 .- 372 .) = 


-f 


y/2m^ / .iL 


\/? 


(A^+A")+4A,x 


y/Xmi ^ p 


\/? 


(A4+Af)-f4A,|| 


(A.5) 


The function ji^iq^) = 1 — 4m^/M tends to 1 for mg <C \/M- 
This condition is well fulfilled for £ = e and q^ ^ I GeV^, 
provided that tensor and scalar Wilson coefficients do not 
receive additional suppression factors. For £ = /r, the value 
of q^ should not be too low, whereas in the case £ = T, 
these observables are not anymore dominated by tensor and 
scalar contributions alone, and the full lepton-mass depen¬ 
dence has to be taken into account. Finally, we note that the 
second part of (7ie -l-72c) in (A.4), 


2mi 


\/M 


At —Ap 


+ 


+ Pl\As\^=4N^XAl^l3i 

2\ 




tiih+ms 


'^P ^P , Cfi' \ 

, I FV®10~ ®10) 

lllh -f Ms q 


(A.6) 


:^{Bs ££) 


1 -t-ys’S^Ar 

1-y? 


£§{Bs -)■ ££). 


(A.8) 


Here y* = Ars/{2rs) with the numerical value given in [28]. 

is the CP asymmetry due to nonvanishing width dif¬ 
ference, which is siX/^r = 1 in the SM, but in general can be 
sX^r G [—1, !]■ Since siX^p can depart from it’s SM value 
in scenarios of new physics considered in this work, we 
take this effect into account in our numerical analysis, al¬ 
though it is suppressed by small y^. The latest SM predic¬ 
tion 'W{Bs ^ Ab) = (3.65 ±0.23) • IQ-*^ [41] includes NLO 
electroweak [42] and NNLO QCD coiTections [43]. 

Finally we discuss the possibility of suitable normal¬ 
izations of 76c, (7ii — 372.) and (7ic ±72c) at high q^ that 
would provide optimized observables. For this purpose we 
use form-factor relations at leading order in 1 /m^ and ne¬ 
glect terms suppressed by mg / With the notation and 

expressions derived in [5], 


7i.-372. = -pf(/i±/jf), 


7lc ± 72e — 3 Pi^/o ± 


3N^X 


{mb-fnisY ° 


(A.9) 


Here /x,||,o and Aq denote B ^ K* form factors, whereas the 
Pj^ depends on vector couplings and pf ± 

Concerning (71^ — 372.), there are no appropriate nor¬ 
malizations, unless the chirality-flipped jq, = 0 be- 






































12 


cause then = Pi holds. There are three potential nor¬ 
malizations 

= (3Pi^+pf )/i + (3Pi +Pi)f^, 

5*.=(pr-pr)/i+(pr-pr)/ii. 

Jis + 2J2s = 3{p^f^+ p^ f«) 


where the last one depends only on vector couplings; i.e., 
is free of tensor and scalar ones. In Jis and J 2 s, tensor cou¬ 
plings contribute either cumulatively or destructively to vec¬ 
tor couplings in p^. 

A similar situation arises for (71,. +J 2 c )7 where no appro¬ 
priate normalization exists, unless scalar couplings vanish. 
In this special case, both Jic and J 2 c = —3/2{p^ — pl)fQ 
depend only on /o and can be used as normalization. In the 
general case, J 2 c still depends only on tensor couplings and 
was used at low [9]. Finally we note that 


Jic-Jic = 3pi /o 
X 


3N^?i 


{mh + msY ** 


(A.ll) 


is free of tensor couplings at high q'^ and would provide ac¬ 
cess to Pj“ provided scalar couplings vanish. 


Appendix B: Theoretical Inputs 

Here we describe the theoretical treatment of observables 
and collect the numerical input for the relevant parameters. 
The software package EOS [12, 40, 44] is used for the cal¬ 
culation of observables in Bg ^ pp and B K^*^£i and 
associated constraints. Both the likelihood and the prior are 
defined entirely within EOS. 

Concerning numerical input, we refer the reader to [28] 
for the compilation of nuisance parameters relevant to this 
work. We adopt the same values for fixed parameters and the 
same priors unless noted otherwise below. Specifically, we 
use identical priors for for common nuisance parameters of 
the CKM quark-mixing matrix, the charm and bottom quark 
masses in the MS scheme, and the parametrization of sub¬ 
leading corrections in 1/mt, as given in [28]. 

Contrary to [28], we do not choose a log-gamma distri¬ 
bution for the asymmetric uncertainties in priors anymore 
but rather a continuous yet asymmetric Gaussian distribu¬ 
tion. This avoids a poor fit because the log-gamma distri¬ 
bution falls off too rapidly in the “short” tail. In a unify¬ 
ing spirit, we now use the continuous rather than discontin¬ 
uous asymmetric Gaussian to approximate asymmetric ex¬ 
perimental intervals in our likelihood. 

The updated prior of the Bg decay constant fs, enters the 
branching ratio of Bg —i' pp. We adopt recent updates of the 


Nf = (2+1) FLAG compilation [45], see Table 5, which av¬ 
erages the results of [46-48]. More recent calculations with 
A/ = (2 -f 1 -f 1) [49] and Nf =2 [50] flavors are consistent 
with these averages. 

The tensor and scalar amplitudes in B factorize 

naively, i.e. they depend only on scalar and tensor B K^*'> 
form factors. These amplitudes are implemented in EOS for 
B Kit and B —>• K*ll as given in [4] and [5], respectively 
and we refrain from using form-factor relations at low and 
high (Y for the tensor and scalar B ^ K and B ^ K* form 
factors. 

Therefore, we need additional nuisance parameters for 
the complete set of B ^ K form factors /+,r.o- As a conse¬ 
quence of using the z parametrization [51] and the kinematic 
relation /o(0) = /+(0), five nuisance parameters (listed in 
Table 5) are needed. As prior information, we average the 
LCSR results [51] and [52] (see Table 5 and cf [28]) sup¬ 
plement them with lattice determinations [13]. The lattice 
results are given in a slightly different parametrization ne¬ 
cessitating a conversion to the parametrization used here. 
For this purpose, we generate the form factors from the 
parametrization of [13], including full correlations, at three 
values of q^ = 17, 20, 23 GeV^. Subsequently, these “data” 
are included in the prior by means of a multivariate Gaus¬ 
sian whose mean and covariance is given in Table 6. The q^ 
values and number of points are chosen such that the cor¬ 
relation of neighboring points is small enough to keep the 
covariance nonsingular. 

We change the parametrization of B ^ K* form factors 
w.r.t. our previous work [28] slightly and now use the sim¬ 
plified series expansion (SSE) [53] 

W)-Y0)]'‘ , (B.l) 

with three parameters {k — 0,1,2) P^r form factor i = 
V, Ao,Ai,Ti,T 23 and two (k = 1,2) for i = Ai2,72. The pa¬ 
rameters Uq correspond to Fi(q^ = 0) and the kinematic con¬ 
straints 

Ai 2(0) = ^kinAo(O), ^[(O) = 72(0), (B.2) 

are used to eliminate where the kinematic factor 

{^MbMk*) depends on the B- and K*- 
meson masses. The pole masses Mr i are set to the values 
given in [23]. 

This parametrization allows us to consistently combine 
available results of form factors from different nonpertur- 
bative approaches, namely LCSR’s at large recoil and lat¬ 
tice at low recoil and to simultaneously implement the kine¬ 
matic constraints. We determine the parameters af in a com¬ 
bined fit to form-factor predictions of V,Ao,i,2,7i 2,3 from 
the LCSR [23] and of y,Ao,Ai,Ai2,71,72,723 from the lat¬ 
tice [24, 25] calculations, including their correlations. We 
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[GeV^] 

17 

20 

23 


0.62 ±0.03 

0.72 ±0.03 

0.87 ±0.04 

f+{q^) 

1.13±0.05 

1.63 ±0.07 

2.68±0.13 


1.02 ±0.06 

1.47 ±0.08 

2.42±0.18 


q^ [GeV2] 

17 

20 

23 

17 

20 

23 

17 

20 

23 

/o(17) 

1.00 

0.93 

0.73 

0.58 

0.48 

0.23 

0.13 

0.010 

0.048 

/o(20) 

- 

1.00 

0.87 

0.41 

0.45 

0.24 

0.059 

0.093 

0.075 

/o(23) 

- 

- 

1.00 

0.28 

0.34 

0.32 

0.010 

0.046 

0.058 

/r(17) 

- 

- 

- 

1.00 

0.78 

0.30 

0.36 

0.21 

0.051 

/+(20) 

- 

- 

- 

- 

1.00 

0.71 

0.22 

0.26 

0.20 

/+(23) 

- 

- 

- 

- 

- 

1.00 

-0.025 

0.14 

0.22 

/r(17) 

- 

- 

- 

- 

- 

- 

1.00 

0.79 

0.46 

fr {20) 

- 

- 

- 

- 

- 

- 

- 

1.00 

0.85 

It {22) 

- 

- 

- 

- 

- 

- 

- 

- 

1.00 


Table 6 Mean values, standard deviations (top) and correlation coefficients (bottom) of lattice points [13] for the 6 —> A" form factors /o,+.r(?^)- 


verified that the constraint 

^12(^Lx) = =^kinAl(^Lx) (B-3) 

at the kinematic endpoint = (Mb—Mk *is satisfied to 
high accuracy by the above constraints and thus need not be 
imposed explicitly. Uninformative flat priors are chosen for 
all a[ with ranges 

c4g[ 0,1], aiG[-2,2], a'e[-3,3]. (B.4) 

The results of the LCSR form-factor predictions have 
been provided to us directly by the authors of [23] at 
= (0.1,4.1,8.1,12.1)GeV2 for VMAiAnJiJiJi 
including the 28 x 28 covariance matrix. Similar results 
could have been obtained by “drawing” form factors from 
the correlated parameters given in [23] and ancillary files 
just as for B ^ K lattice form factors. 

For the B ^ K* lattice form factors this approach is not 
good enough as we were not able to select more than two 
values without obtaining a singular covariance. In order to 
fully the exploit the available information, we then contacted 
the authors of [25] and obtained the original values of the 
form factors (including correlation) at various values in the 
interval q^ ^ [11-9, 17.8] to which Horgan et al. fit the SSE. 
The covariance has a block-diagonal structure with a 48 x 48 
block for V,Ao,Ai,Ai 2 and a 36x36 block for 71,72,723. 
Having the “raw” information on form-factor values is much 
more reliable and future proof as there are no issues with 
artificial correlation and we could one day decide to use yet 
another form-factor parametrization and fit it easily to these 
data points. 

We have compared the SSE fit (B.l) with two versus 
three parameters and found that in the former case lattice 


form factors influence the fit such that form factors tend to 
be higher than LCSR predictions at low q^ leading to a poor 
fit. Hence we prefer the three-parameter setup as it provides 
the flexibility needed to accommodate LCSR and lattice re¬ 
sults. Means and standard deviations of that three-parameter 
fit are given in Table 5, we omit the correlations for the sake 
of brevity but are happy to provide them. 

Appendix C: Monte Carlo sampling 

The marginalization of the posterior is performed with 
the package pypmc [54], which incorporates the algo¬ 
rithm presented in [55, 56] and in addition an implemen¬ 
tation of the variational Bayes algorithm. In every analy¬ 
sis we first run multiple adaptive Markov chains (MCMC) 
in parallel through pypmc. If necessary, chains are seeded 
at the SM point to exclude solutions in which multiple 
nuisance parameters—mostly for hadronic corrections— 
simultaneously deviate strongly from prior expectations. 

In total, there are 19 parameters aj to describe B ^ K* 
form factors and most of them are strongly correlated. But it 
is well known that strong correlation leads to poor sampling 
as it can cause the random-walk Markov chains to spend an 
excessive amount of time in regions of low probability and 
thus produce spurious peaks. To mitigate this issue, we per¬ 
form a fit to form-factor constraints without any experimen¬ 
tal data and use the resulting covariance matrix to transform 
parameters such that the new parameters are uncorrelated. 

In all but three cases, the Markov chains then give reli¬ 
able results. But when we analyze scenarios with '^S,P 7^ ® 
and all experimental constraints, strong correlations appear 
again. As a final solution, we then use importance sampling 
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Quantity 

Prior 

Reference 

fi, decay constant 

/fi. 

(227.7 ±4.5) MeV 

[45] 

B —> A" form factors 

/+(0)=/o(0) 

0.34 ±0.05 

[51,52] 

/r(0) 

0.38 ±0.06 

[51,52] 


—4 

[51] 

bi 

-2 

[51] 


o 9+1.0 
^•^-2.0 

[51] 

B ^ K* fomi factors 

Uq 

0 3 ^+o .02 

[23, 25] 


1 9 1 +0.08 
^•^^-0.08 

[23, 25] 


0 77+0.45 
' '-0.48 

[23, 25] 

«o‘ 


[23, 25] 


0-0lt°om 

[23, 25] 

«2^‘ 

0.045+°;|2 

[23, 25] 


0 53+0-10 
^•‘^'^-0.14 

[23, 25] 


^•■^^-0.39 

[23, 25] 

< 

0 35+0.02 
^••5J_o.03 

[23, 25] 

«r 

_ 1 1 q+0.08 
^•*^-0.08 

[23, 25] 


1 55+0.33 

*•‘^'^-0.33 

[23, 25] 


0 31+0.02 

^•■^*-0.02 

[23, 25] 


i-U/_o.o8 

[23, 25] 

al' 

*•^^-0.30 

[23, 25] 

al^ 

n 99+0.08 
^••5J_0.10 

[23, 25] 

al^ 

^•^^-0.24 

[23, 25] 

«o“ 

^•°'-0.05 

[23, 25] 


0 03+0.28 
^■^”^-0.22 

[23, 25] 

al^^ 

—0 15+0-80 
‘4-i2)_0 7i 

[23, 25] 


Table 5 Prior distributions of the nuisance parameters for hadronic 
quantities. 


with the initial proposal function determined by a fit of a 
Gaussian mixture to the MCMC samples within the vari¬ 
ational Bayes approximation [57]. As the posterior is uni- 
modal and closely resembles a Gaussian, only a few Gaus¬ 
sian components are needed; i.e., 3 components proved opti¬ 
mal by the variational approximation to the model evidence. 

In the most challenging run with 62 parameters, we ob¬ 
tain a relative effective sample size of only 0.038%. We want 
to have enough independent samples N such that the 68% re¬ 
gion is determined with a relative precision of about 1%. As 
a rule of thumb, we consider the “relative error of the error” 
given by 1 /'/TN [38, ch. 37]. In the 62D case, we compute 
a total of 1.1 • 10® importance samples, update the proposal 
after every 10® samples, and combine all samples [57] such 


that N « 3500 and the estimated relative error of the error is 
1.2% and thus good enough for our purposes. 

To create the smooth marginal plots in Figures 1-3, we 
apply kernel density estimation for both MCMC and impor¬ 
tance samples using the fast figtree library [58]. In the 
latter case, we additionally crop 500 outliers. 

The prior or posterior predictive distribution of an ob¬ 
servable X within a model M in which X — f {9) is a definite 
function of the parameters 9 and given as 

P{X\M) = J d9P{X\9,M)P{9\M) 

= Jd9 5{X-fi9))Pi9\M). (C.l) 

We estimate P{X\M) by computing /(0,) for every sample 
9i ~ P{9\M), then smooth as above. 
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